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ABSTRACT 

We report the discovery of HAT-P-17b,c, a multi-planet system with an inner transiting planet 
in a short-period, eccentric orbit and an outer planet in a 4.8 yr, nearly circular orbit. The inner 
planet, HAT-P-17b, transits the bright V = 10.54 early K dwarf star GSC 2717-00417, with an 
orbital period P = 10.338523 ± 0.000009 d, orbital eccentricity e = 0.346 ± 0.007, transit epoch 
T c = 2454801. 16945±0.00020 (BJE0), and transit duration 0.1691 ± 0.0009d. HAT-P-17b has a mass 
of 0.530 ±0.018 Mj and radius of 1.010 ±0.029 Rj yielding a mean density of 0.64 ± 0.05 g cm" 3 . This 
planet has a relatively low equilibrium temperature in the range 780-927 K, making it an attractive 
target for follow-up spectroscopic studies. The outer planet, HAT-P-17c, has a significantly longer 
orbital period P 2 = 17981J9 d and 

a minimum mass 7712 sin 12 = L4l ! 4 M3. The orbital inclination of 
HAT-P-17c is unknown as transits have not been observed and may not be present. The host star has 
a mass of 0.86 ± 0.04 M©, radius of 0.84 ± 0.02 i?©, effective temperature 5246 ± 80 K, and metallicity 
[Fe/H] = 0.00 ±0.08. HAT-P-17 is the second multi-planet system detected from ground-based transit 
surveys. 

Subject headings: planetary systems — stars: individual (HAT-P-17, GSC 2717-00417) techniques: 
spectroscopic, photometric 



1. INTRODUCTION 

With nearly 100 confirmed transiting extrasolar plan- 
ets (TEPs) known, many studies of planetary properties 
now focus on the statistical distributions of and corre- 
lations between planetary parameters. Individual TEPs 
still remain extraordinarily valuable, particularly if they 
have properties that exemplify an important subgroup of 
planets and orbit stars that are bright enough for mean- 
ingful follow-up observations. Such iconic well-studied 
planets include HP 209458b dCharbonneau et al.l l2000t 
IHenrv et all [2000h. HP 189733 (IBouchv et alJ 120051 
GJ 43 6b (IGillon et , allllM IButler et al.ll2004D. HAT-P - 
13b,c (Ba kos et al.ll2010fl. WASP-12b (IHebb et al.| [2009L 
and GJ 1214b (|Charbonneau~et~aT1l2009D . The HAT-P- 
17 system has at least two unusual properties compared 
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to the ensemble of known TEPs and may serve as an 
exemplar for planets with these properties. The atmo- 
sphere of HAT-P-17b is relatively cool for a TEP and 
HAT-P-17c is one of only two long-period planets found 
to orbit a TEP host. 

The distribution of TEPs discovered by ground-based 
transit surveys is biased toward large planets orbiting 
faint early type stars in short period orbits. Each of 
these biases stems from the observational selection ef- 
fects of the surveys that have detected the majority of 
TEPs: the deep transits of large planets are easier to 
detect; early type stars dominate magnitude limited sur- 
veys; and short-period orbits have higher a prioiri transit 
probabilities and significantly larger probabilities of de- 
tection in a ground-based survey limited to one or two 
observing sites. The overabundance of detected short 
period planets has skewed our perception of the atmo- 
spheric properties of extrasolar jovian planets. Plan- 
ets in F ~ 3d orbits (a ~ 0.04) experience intense 
interactions with the radiation and tides of their host 
stars. The atmospheres of many of these planets are 
infl ated beyond the rad ii predicted by theoretical mod- 
els (jFortnev et al.| [2007). To understand cooler planets, 
which we know from radial velocity (RV) surveys repre - 
sent the vast majority of gas giants (Wri ght et al.l l2009). 
we must study planets orbiting progressively further from 
their host stars. 

With such a scarcity of cool planets, HAT-P-17b is 
a valuable probe of the planetary mass-radius relation- 
ship and additional properties through follow-up ob- 
servations. Together, the relatively long orbital pe- 
riod and later spectral type of the host star yield an 
incident stellar flux that is 1-2 orders of magnitude 
lower than the flux received by most detected TEPs 
(|Kovacs et al.l I2010D . The host star is also relatively 
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bright, V — 10.54, making follow-up atmo spheric studies 
conce ivable. While the Kepler mission ([Borucki et al.l 
I2010T) has been extraordinarily successful in the detec- 
tion of TEPs, the vast majority of its discovered planets 
orbit f aint stars; only 1.5% of the ~ 10 5 stars being sur- 
veyed (Batalh a et al.lr2~01Q[ ) are brighter than 11th mag- 
nitude (Kepler magnitude). We predict that HAT-P-17b 
will be among the small number of well-studied cooler 
(T cff < 1000 K) TEPs. 

Prior to this announcement, only one TEP discovered 
by a ground-based photometric s urvey is in a confir med 
multi-planet system. HAT-P-13 (jBakos et alJl2010f > has 
a hot Jupiter inner planet and highly eccentric super- 
Jupiter outer planet with an orbital period of 450 d. 
The outer planet, HAT-P-13c, was detected only in RV 
measurements and has not been shown to transit. The 
system reported here, HAT-P-17, is now the second TEP 
discovered by a ground-based transit survey with a con- 
firmed second planet. The relatively low rate of de- 
tected planet multiplicity among TEPs discovered from 
by ground may be skewed by the lack of long-term RV 
and photometric monitoring for most TEP host stars. 
Measuring the rate of planet multiplicity among stars 
hosting a hot Jupiter will probe the dynamical histories 
and migrations mechanisms of hot Jupiters (see, e.g., 
IWu et all 120071 ) . Multi-planet systems ar e significantly 
more c ommon among RV-detected systems: iWright et al.l 
(2009) find that 28% of known planet host stars are 
multi-planet systems. 

Several multi-planet systems have also been discovered 
from space. C orot-7 is thought to host two short-p eriod 
super-Earths (jLeger et. al.ll2009HQueloz et al.ll2009l) . one 
of which transits. The Kepler mission recently an- 
nounced five candidate systems with multiple transit- 
ing planets and is p oised to announce additional systems 
(Stc ffen et al.ll20toh . 

The Hu ngarian- made Aut omated Telescope Network 
(HATNet; IBakos et al"1l2004D survey has been one of the 
main contributors to the discovery of TEPs. In op- 
eration since 2003, it has now covered approximately 
14% of the sky, searching for TEPs around bright stars 
(8 < / < 14.0). HATNet operates six wide-field instru- 
ments: four at the Fred Lawrence Whipple Observatory 
(FLWO) in Arizona, and two on the roof of the hangar 
servicing the Smithsonian Astrophysical Observatory's 
Submillimeter Array, in Hawaii. Since 2006, HATNet 
has announced and published 16 TEPs. In this work we 
report our seventeenth discovery, around the relatively 
bright star previously known as GSC 2717-00417. 

The layout of the paper is as follows. In Section [5] we 
report the detection of the photometric signal and the 
follow-up spectroscopic and photometric observations of 
HAT-P-17. In Section |3] we describe the analysis of the 
data, beginning with the determination of the stellar pa- 
rameters, continuing with a discussion of the methods 
used to rule out nonplanetary, false positive scenarios 
which could mimic the photometric and spectroscopic ob- 
servations, and finishing with a description of our global 
modelling of the photometry and radial velocities. In 
Section |4] we discuss implications of this discovery, com- 
pare our results with recent theoretical models of TEPs, 
and consider possible follow-on observations. 

2. OBSERVATIONS 



2.1. Photometric detection 

The transits of HAT-P-17b were detected with the 
HAT-5 telescope in Arizona, the HAT-8 telescope in 
Hawaii, and with the Wis e-HAT (WHAT) tel escope at 
Wise Observatory in Israel (|Shporer et al.ll2009l ). The re- 
gion around GSC 2717-00417, a field internally labeled as 
247, was observed on a nightly basis between 2005 May 8 
and 2005 October 24, whenever weather conditions per- 
mitted. We gathered 9686 exposures of 5 minutes du- 
ration at a 5.5 minute cadence. Each image contained 
approximately 85,000 stars down to / ~ 14.0. For the 
brightest stars in the field, we achieved a per-image pho- 
tometric precision of 5 mmag. The star is also located in 
the overlapping field 248, which was observed with the 
HAT-6 telescope in Arizona and the WHAT telescope in 
Israel between 2004 June 4 and 2004 November 10, and 
between 2005 July 3 and 2005 July 15. Altogether 4882 
exposures of 5 minutes duration at 5.5 minute cadence 
were gathered for this field. 
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Fig. 1. — Unbinned (top) and binned (bottom) light curves of 
HAT-P-17 including all 14,000 instrumental /-band 5.5 minute ca- 
dence measurements obtained with the HAT-5, HAT-6, and HAT-8 
telescopes of HATNet and with the WHAT telescope (see the text 
for details), and folded with the period P = 10.3385228 days (re- 
sulting from the global fit described in Section |3J. The soli d lin e 
shows the "P1P3" transit model fit to the light curve fSection !3.3l l. 

The calibration of the HATNet and WHAT frames 
was carried out using standard photometric procedures. 
The calibrated images were then subjecte d to star de- 
tectio n and astrometry, as described in IPal fc Bakos 
(2006). Aperture photometry was performed on each 
image at the stellar centroids derived from the Two Mi- 
cron All Sky Survey (2MASS; iSkrutskie et~aLll2006l) cat- 
alog and the individual astrometric solutions. The re- 
sulting light curves were decorrelated (cleaned of trends) 
using the Externa l Parameter Decorrelation (EPD; see 
IBakos et al.l 120101) technique in "constan t" mode and 
the T rend Filtering Algorithm (TEA; see IKovacs et al.l 
2005). The light curves were searched for periodic 
bo x-shaped signals usi ng the Box Least-Squares (BLS; 
see IKovacs et al.l I2002T ) method. We detected a signifi- 
cant signal in the light curve of GSC 2717-00417 (also 
known as 2MASS 21380873+3029193 and TYC 2717- 
417-1; a = 21 h 38 m 08 88s, 6 = +30°29'19.4"; J2000; 
V=10.54; iDroege et al.l [20061 ) . with an apparent depth 
of ~ 12.0 mmag, and a period of P = 10.3385 days 
(see Figure [lj. The drop in brightness had a first-to- 
last-contact duration, relative to the total period, of 
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q = 0.0164 ±0.0001, corresponding to a total duration of 
Pq = 4.058 ± 0.022 hr (see Figure [J). We note that the 
transits were only detected from the observations of field 
247, and were not detected in the observations of field 
248. 

2.2. Reconnaissance Spectroscopy 

As is routine in the HATNet project, all candidates 
are subjected to careful scrutiny before investing valu- 
able time on large telescopes. This includes spectro- 
scopic observations at relatively modest facilities to es- 
tablish whether the transit-like feature in the light curve 
of a candidate might be due to astrophysical phenom- 
ena other than a planet transiting a star. Many of these 
false positives are associated with large radial-velocity 
variations in the star (tens of kms -1 ) that are easily rec- 
ognized. 

One of the tools we have used for this purpose is the 
Harvard-Smithsonian C enter for Astro physics (CfA) Dig- 
ital Speedometer (DS; iLathaml 119921 ) . an echelle spec- 
trograph mounted on the FLWO 1.5 m telescope. This 
instrument delivers high- resolution spectra (A/ A A « 
35,000) over a single order centered on the Mg lb triplet 
(~5187A), with typically low signal-to-noise (S/N) ra- 
tios that are nevertheless sufficient to derive radial veloc- 
ities (RVs) with moderate precisions of 0.5-1.0 km s _1 for 
slowly rotating stars. The same spectra can be used to 
estimate the effective temperature, surface gravity, and 
pr ojected rotati o nal ve locity of the host star, as described 
bv lTorres et alJ (|2002h . With this facility we are able to 
reject many types of false positives, such as F dwarfs or- 
bited by M dwarfs, grazing eclipsing binaries, or triple or 
quadruple star systems. Additional tests are performed 
with other spectroscopic observations described in the 
next section. 

For HAT-P-17 we obtained eight observations with 
the DS between September and November of 2007. 
The velocity measurements showed an rms residual of 
0.60 km s -1 , consistent with no detectable RV variation 
within the precision of the measurements. All spectra 
were single-lined, i.e., there is no evidence for additional 
stars in the system. The atmospheric parameters we 
infer from these observations are the following: effec- 
tive temperature T c ff+ = 5250 ± 125 K, surface gravity 
log g+ = 4.5 ±0.25 (log cgs), and projected rotational ve- 
locity vsini = O.O^q'q kms -1 . The effective temperature 
corresponds to an early K dwarf star. The mean helio- 
centric RV of HAT-P-17 is 7 RV = 20.13 ± 0.21 kms" 1 . 

2.3. High resolution, high S/N spectroscopy 

Given the significant transit detection by HATNet, and 
the encouraging DS results that rule out obvious false 
positives, we proceeded with the follow-up of this can- 
didate by obtaining high-resolution, high-S/N spectra to 
characterize the RV variations, and to refine the deter- 
minatio n of the stellar parameters. For this we used 
HIRES (|Vogt et al.lll994D on the Keck I telescope located 
on Mauna Kea, Hawaii, between 2007 October and 2010 
April. The width of the spectrometer slit was 0'.'86, re- 
sulting in a resolving power of A/AA « 55,000, with a 
wavelength coverage of -3800-8000 A. 

We obtained 32 HIRES exposures with an iodine cell 
mounted directly in front of the spectrometer entrance 



slit. The dense set of molecular absorption lines im- 
printed on the stellar spectra provide a robust wave- 
length fiducial against which Doppler shifts are mea- 
sured, as well as strong constraints on the shape of the 
spectrome ter instrumental profile at the time of eac h 6b- 
servation ([Marcv fc ButleJl992tl"Valenti et aDl995[ ). An 
additional exposure was taken without the iodine cell, for 
use as a template in the reductions. Relative RVs in the 
so lar system ba r ycent ric frame were derived as described 
bv lButler et al.l (|l996| ) , incorporating full modeling of the 
spatial and temporal variations of the instrumental pro- 
file. These measurements have typical uncertainties of 
1.5-2.0 ms _1 for spectra with per-pixel signal-to-noise 
ratios of 100-150. HIRES measurements of late G and 
early K dwarf stars have achieved long term stability of 
1.5-2.0 ms _1 on standard stars, inc luding noise from sys - 
tematic and astrophysical sources (jHoward et alJ[20Ta ). 
The RV measurements and their uncertainties are listed 
in Table [TJ The period- folded data, along with our best 
fit described below in Section [3l are displayed in Figure[2j 

In the same figure we also show the relative S index, 
which is a measure of the chromospheric activity of the 
star derived from the flux in the cores of the Ca II H and 
K lines. This i ndex was computed following the pre scrip- 
tion given by iVaughan. Preston fc Wilson! (|1978l ) , after 
matching each spectrum to a reference spectrum using 
a transformation that includes a wavelength shift and a 
flux scaling that is a polynomial as a function of wave- 
length. The transformation was determined on regions 
of the spectra that are not used in Note that our rel- 
ative S index has not been calibr ated to the scale of 
IVaughan. Preston fe Wilson! ()1978[ ). We do not detect 
any significant variation of the index correlated with the 
orbital phase of either planet; such a correlation might 
have indicated that the RV variations could be due to 
stellar activity, casting doubt on the planetary nature of 
the candidates. 

In addition, we computed an Shk index calibrated 
to the Mt. Wilson scale, permitting comparisons 
with calibrated activity measureme nts of other stars 
(jKnutson. Howard, fc Isaacson! 12010). We find S'hk = 
0.162 (median of all HIRES measurements) and \ogR^ K 
= —5.039 (median). These measurements empl oy the 
techniques described in llsaacson fc Fischer! (|2010l ). cali- 
brated on 1500 stars observed by the California Planet 
Survey (CPS). We used B—V = 0.83 estimated from T cff 
us ing the linear tran s forma tion between those quantities 
in lValenti fc Fischer! (f2005h . 

2.4. Photometric follow-up observations 

In order to permit a more accurate modeling of the 
light curve, we conducted additional photometric obser- 
vations using a variety of facilities, including: the Ke- 
plerCam CCD camera on the FLWO 1.2 m telescope, the 
0.6 m Schmidt telescope of Konkoly Observatory at the 
Piszkesteto Mountain Station, and the 0.46 m and 1.0 m 
telescopes at Wise Observatory. We observed three tran- 
sit events of HAT-P-17 with the FLWO 1.2 m telescope 
on the nights of 2007 December 14, 2008 October 19, 
and 2009 October 16, while the event on 2008 September 
8 was observed simultaneously at Konkoly Observatory 
and with the two telescopes at Wise Observatory (Fig- 
ure [3]) . These observations are summarized in Table [2] 

The reduction of these images, including basic cali- 
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TABLE 1 

Relative radial velocities, bisector spans, and activity 
index measurements of hat-p-17. 
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Note. — Note that for the iodine-free template exposures we 
do not measure the RV but do measure the BS and S index. Such 
template exposures can be distinguished by the missing RV value. 
a The zero-point of these velocities is arb itrary. An overall offset 
7 ro l fitted to these velocities in Section l3.3l has not been subtracted. 
b Internal errors ex cludi ng the component of astrophysical jitter 
considered in Section 13.31 

c Relative chromospheric activ i ty ind ex, not calibrated to the scale 
of lVaughan. Preston fc Wilsonl l|l978T l. 



TABLE 2 

Summary of photometric follow-up observations 
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Fig. 2. — Top: Keck/HIRES RV measurements for HAT-P-17 shown as a function of BJD, along with our best-fit 2-planet model (see 
TablefSJ. The center-of-mass velocity has been subtracted. Second from top: Same as top panel except the RV model of the inner planet has 
been subtracted from the data and the model, revealing the orbit of the outer planet. The rms variation of the residuals to the two-planet 
model is 3.07ms -1 , requiring a jitter of 2.3ms -1 added in quadrature to the individual errors to yield a reduced \ 2 of 1.0. The error-bars 
in this panel have been inflated accordingly. Third row: RV measurements phased to the orbital periods of the inner planet (left) and the 
outer planet (right). In each plot the orbit of the other planet has been removed. Fourth row: Bisector spans (BS), with the mean value 
subtracted, phased a t the period of the inner planet (left) and the outer planet (right). The measurement from the template spectrum is 
included (see Section l3.2H . Fifth row: Relative chromospheric activity index S measured from the Keck spectra, phased at the period of the 
inner planet (left) and the outer planet (right). Note the different vertical scales of the panels. Observations shown twice are represented 
with open symbols. 
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TABLE 3 

High-precision differential photometry of HAT-P-17 



-0.1 -0.05 0.05 0.1 
Time from transit center (days) 

Fig. 3. — Unbinned instrumental 2-band, /-band, i?-band, 
and i-band transit light curves, acquired with KeplerCam at the 
FLWO 1.2 m telescope, with the Schmidt telescope at Konkoly Ob- 
servatory, and with the 0.46 m and lm telescopes at Wise Obser- 
vatory. The li ght curves have been EPD- and TFA-processed, as 
described in S 13.31 The dates and facilities used to observe the 
events are indicated. Curves after the first are displaced verti- 
cally for c larit y. Our best fit from the global modeling described 
in Section 13.31 is shown by the solid lines. Residuals from the fits 
are displayed at the bottom, in the same order as the top curves. 
The error bars represent photon and background shot noise, plus 
readout noise. 

bration, astrometry, and aperture photomet ry, was per- 
formed as described by iBakos et al.l (|2010T ) . We per- 
formed EPD and TFA to remove trends simultane- 
ously with the ligh t curve modeling (see Section |31 and 
IBakos et all ([2010) for details). The final time scries are 
shown in the top portion of Figure [3l along with our 
best-fit transit light curve model described below; the 
individual measurements are reported in Table |3] 

3. ANALYSIS 

3.1. Properties of the parent star 

Fundamental parameters of the host star HAT-P-17 
such as the mass (M*) and radius (-R*), which are needed 
to infer the planetary properties, depend strongly on 
other stellar quantities that can be derived spectroscop- 
ically. For this we have relied on the HIRES tem- 
plate spectrum, and the anal ysis package known as Spec - 
troscopy Made Easy (SME; iValenti fc Piskunovl fl 996), 
along with the atomic line database of l Valenti fc FischerJ 
( 2005) . SME yielded the following initial values and un- 
certainties (which we have conservatively increased to in- 
clude our estimates of the systematic errors): effective 
temperature T e ff* = 5455 ± 80 K, stellar surface gravity 
logs* = 4.60 ± 0.10 (cgs), metallicity [Fe/H] = +0.13 ± 
0.08 dex, and projected rotational velocity vsini = 
0.5 ± 0.5 kms -1 . We adopt the single-sided uncertainty 



BJD 

(2,400,000+) 


Mag a 


"Mag 


Mag(orig) b 


Filter 


54449.55537 


-0.00357 


0.00095 


9.35242 


z 


54449.55568 


0.00062 


0.00094 


9.35614 


z 


54449.55603 


-0.00175 


0.00095 


9.35430 


z 


54449.55633 


0.00203 


0.00095 


9.35788 


z 


54449.55668 


-0.00602 


0.00093 


9.35026 


z 


54449.55700 


0.00173 


0.00095 


9.35703 


z 


54449.55732 


-0.00628 


0.00093 


9.34882 


z 


54449.55768 


-0.00183 


0.00094 


9.35300 


z 


54449.55799 


-0.00084 


0.00094 


9.35409 


z 


54449.55835 


0.00139 


0.00093 


9.35711 


z 



Note. — This table is available in a machine- re ad able form in the 
online journal. A portion is shown here for guidance regarding its form 
and content. 

a The out-of-transit level has been subtracted. These magnitudes have 
been subjected to the EPD and TFA procedures, carried out simulta- 
neously with the transit fit. 

k Raw magnitude values without application of the EPD and TFA 
procedures. 



of ±0.5 km s on ^sinz from IValenti fc Fischer! ([2005) 
based on their SME analysis of nearly 2000 stars. For 
this star and others with low usini, the true error dis- 
tribution excludes unphysical values (usini < Okms -1 ) 
and is likely asymmetric. 

In principle the effective temperature and metallicity, 
along with the surface gravity taken as a luminosity in- 
dicator, could be used as constraints to infer the stellar 
mass and radius by comparison with stellar evolution 
models. However, the effect of logg* on the spectral 
line shapes is subtle, and as a result it is typically dif- 
ficult to determine accurately, so that in practice it is 
a poor luminosity indicator. For planetary transits, a 
stronger constraint is often provided by a/R*, the nor- 
malized semimajor axis, which is closely related to p*, the 
mean stellar density. The quantity a/R* can be derived 
direc tly from the transit light curves (see ISozzetti et all 
I2007L and also Section l3~3| . This, in turn, improves our 
determination of the spectroscopic parameters by sup- 
plying an indirect constraint on the weakly determined 
spectroscopic value of log<?* that removes degeneracies. 
We take this approach here, as described below. The 
validity of our assumption, namely that the adequate 
physical model describing our data is a planetary transit 
(as opposed to a blend) , is shown later in Section 13.21 

Our initial values of T ff+, logg*, and [Fe/H] were 
used to determine auxiliary quantities needed in the 
global modeling of the follow-up photometry and radial 
velocities (specifically, the limb-darkening coefficients). 
This modeling, the details of which are described in 
Section |3"U1 uses a Monte Carlo approach to deliver the 
probability dis t ribution of a/R+ and other fitted vari- 
ables. See EH (|2009bft for further details. When com- 
bining a/ R+ (a luminosity proxy) with assumed Gaussian 
distributions for T e fT* and [Fe/H] from SME, we compare 
with stellar evolution models to estimate the probabil- 
ity distributions of additional inferred stellar parameters, 
including logg*. Here we use the stellar evol ution cal- 
culati ons from the Yonsei-Yale (YY) series by lYi et al.1 
pOOl . The comparison with the model isochrones was 
carried out for each of 20,000 Monte Carlo trial sets (see 
Section |373| . Parameter combinations corresponding to 
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unphysical locations in the H-R diagram (1.5% of the tri- 
als) were ignored, and replaced with another randomly 
drawn parameter set. The result and error estimate for 
the surface gravity, log g* = 4.53 ±0.02, is different from 
the result of our initial SME analysis, which is not sur- 
prising in view of the strong correlations among Teg*, 
[Fe/H], and log that are often present in spectroscopic 
determinations. Therefore, we carried out a second iter- 
ation in which we adopted this value of log g+ and held it 
fixed in a new SME analysis (coupled with a new global 
modelling of the RV and light curves), adjusting only 
T off *, [Fe/H], and vsini. This gave T cff * = 5246 ± 80 K, 
[Fe/H] = 0.00 ± 0.08, and vsini = 0.3 ± 0.5 km s" 1 , in 
which the conservative uncertainties for the first two have 
been increased by a factor of two over their formal val- 
ues, as before. A further iteration did not change log g* 
significantly, so we adopted the values stated above as 
the final atmospheric properties of the star. They are 
collected in Table IU together with the adopted values 
for the macroturbulent and microturbulent velocities. 

With the adopted spectroscopic parameters the model 
isochrones yield the stellar mass and radius, M* = 
O.857±O.O39M and = 0.837± 0.021 R , along with 
other properties listed at the bottom of Table [4] HAT- 
P-17 is a early K dwarf star with an estimated age of 
7.8 ± 3.3 Gyr, according to these models. The inferred 
location of the star in a diagram of o/i?* versus T e ff*, 
analogous to the classical H-R diagram, is shown in Fig- 
ure |U The stellar properties and their la and 2a con- 
fidence e l lipsoi ds are displayed against the backdrop of 
lYi et al.l (|2001| ) isochrones for the measured metallicity 
of [Fe/H] = 0.00, and a range of ages. For comparison, 
the location implied by the initial SME results is also 
shown (triangle). 




Effective temperature [K] 

Fig. 4. — Model isochrones from Yi ct al. (2001) for the measured 
metallicity of HAT-P-17, [Fe/H]= 0.00, and ages between 1.0 and 
13.0 Gyr with a step-size of 1.0 Gyr (left to right). The adopted 
values of and a/i?* are shown together with their lrr and 2cr 

confidence ellipsoids. The initial values of T^ff* and a/R+ from the 
first SME and light curve analyses are represented with a triangle. 

The stellar evolution modelling provides color indices 
that may be compared against the measured values as 
a sanity check. The best available measurements are 
the ne ar-infrared magnitud es from the 2MASS Cata- 
logue (jSkrutskie et all [2006h . J 2M ass = 9.017 ± 0.022, 
-H2MASS = 8.619 ± 0.029 and K 2MAS s = 8.544 ± 0.025; 
which we have converted to the photometric system of 
the models (ESO system) using the transformations by 



TABLE 4 
Stellar parameters for HAT-P-17 



Parameter Value Source 



Spectroscopic properties 

T cB * (K) 5246 ± 80 SME a 

[Fe/H] 0.00 ± 0.08 SME 

vsini (kms" 1 ) . . . 0.3 ± 0.5 SME 

u mac (kms" 1 ) . . . . 3.21 SME 

v mic (kms' 1 ).... 0.85 SME 

Trv (kms" 1 ) 20.13 ±0.21 DS 

Photometric properties 

V (mag) 10.54 TASS 

V- I a (mag) 0.901 ±0.10 TASS 

J (mag) 9.017 ±0.022 2MASS 

H (mag) 8.619 ± 0.029 2MASS 

K a (mag) 8.544 ± 0.025 2MASS 

Derived properties 

M* (M@) 0.857 ±0.039 YY+a/i?*+SME b 

R* (R e ) 0.837 ±0.021 YY+a/R*+SME 

logg* (cgs) 4.53 ±0.02 YY+a/R*+SME 

L* (£0 ) 0.48 ±0.04 YY+a/i?*+SME 

M v (mag) 5.75 ±0.12 YY+a/i?*+SME 

M K (mag,ESO) . . 3.79 ± 0.07 YY+a/fi*+SME 

Age (Gyr) 7.8 ± 3.3 YY+a/i?*+SME 

Distance (pc) 90 ± 3 YY+a/R^+SME 



a SME — "Spectroscopy Made Easy" package for the analysis of 
high-resolution spectra (Valcnti & Piskunov 1996). These param- 
eters rely primarily on SME, but have a small dependence also on 
the iterative analysis incorporating the isochrone search and global 
modelling of the data, as described in the text. 

b YY+o/fl,+SME = Based on the YY isochrones lYi et all20uTT ), 
a/-R* as a luminosity indicator, and the SME results. 

ICarpenterl (|2001l ) . The resulting measured color index is 
J-K = 0.504 ±0.036. This is within la of the predicted 
value from the isochrones of J — K = 0.53 ± 0.02. The 
distance to the object may be computed from the abso- 
lute K magnitude from the models (Mr = 3.79 ± 0.07) 
and the 2MASS K s magnitude, which has the advantage 
of being less affected by extinction than optical magni- 
tudes. The result is 90 ± 3 pc, where the uncertainty ex- 
cludes possible systematics in the model isochrones that 
are difficult to quantify. 

3.2. Spectral line-bisector analysis 

Our initial spectroscopic analyses discussed in 
Section l2~2l and Section l2~3l rule out the most obvious as- 
trophysical false positive scenarios. However, more sub- 
tle phenomena such as blends (contamination by an un- 
resolved eclipsing binary, whether in the background or 
associated with the target) can still mimic both the pho- 
tometric an d spectroscop i c sign atures we see. 

Following iTorres et al.1 (|2007l ). we explored the possi- 
bility that the measured radial velocities are not the re- 
sults of a planet in Keplerian motion, but are instead 
caused by distortions in the spectral line profiles due to 
contamination from a nearby unresolved eclipsing binary. 
A bisector analysis b ased on the Keck sp ectra was done 
as described in §5 of iBakos et al.l (|2007l ). We detect no 
variation in excess of the measurement errors in the bi- 
sector spans (see Figure [2|) . The correlation between the 
radial velocities and the bisector variations is insignifi- 
cant. Therefore, we conclude that the velocity variations 
are real, and that the star is orbited by a close-in giant 
planet. 
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3.3. Global modelling of the data 

This section describes the procedure we followed to 
model the HATNet photometry, the follow-up photome- 
try, and the radial velocities simultaneously. Our model 
for the follow-up light curves used analytic formulae 
based on lMandel fc Ago! ()2002l ) for the eclipse of a star 
by a planet, with limb darkening being prescribed by a 
quadratic law. The limb darkening coefficients for the 
/, R and Sl oan z and i b ands were interpolated from 
the tables bv lClaretl (|2004) for the spect roscopic param- 
eters of the star as determined from the SME analysis 
f Section 13. ip . The transit shape was parametrized by 
the normalized planetary radius p = R p / R+, the square 
of the impact parameter b 2 , and the reciprocal of the half 
duration of the transit C/-R*. We chose these parame- 
ters because of their simple geometric meani ngs and their 
negli gible correlations with each other (see iBakos et al.l 
2010D . The relation between (,/R* and the quantity 
a/R*, used in Section |3~TI is given by 

a/R+ = P/2tt(C/-R*)V1 - 6Vl - e 2 /(l + esinw) (1) 

(see, e.g., iTinglev fc Sa ckctt 20Q3). Our model for the 
HATNet data was the simplified "P1P3" version of the 
iMandel &: Agoll (|2002f) analytic functions (an expansion 
in terms of Legendre polyno mials), for the reasons de- 
scribed in IBakos et al.l (|2010T ) . 

Initial modelling of the RV observations showed devi- 
ations from a Keplerian fit highly suggestive of a sec- 
ond body in the system with a much longer period 
than the transiting planet. Thus, in our global mod- 
elling, the RV curve was parametrized by the combina- 
tion of an eccentric Keplerian orbit for the inner planet 
with semi-amplitude K, and Lagrangian orbital elements 
(k,h) = e x (cosw, sinw), plus an eccentric Keplerian 
orbit for the outer object with K~2, ki and /12, and a 
systemic RV zero-point 7 (see also Ba kos et all 120101 ). 
Throughout this paper the subscripts "1" and "2" re- 
fer to HAT-P-17b and HAT-P-17c, respectively. If the 
subscript is omitted, we refer to HAT-P-17b. 

We assumed a strict periodicity in the individual tran- 
sit times. We assigned the transit number N tr = to the 
last complete follow-up light curve gathered on 2009 Oct 
16. The adjustable parameters in the fit that determine 
the ephemeris were chosen to be the time of the first 
transit center observed with HATNet, T Ci _65, and that 
of the last transit center observed with the FLWO 1.2 m 
telescope, T Cy0 . We used these as opposed to period and 
reference epoch in or der to minimize correlations be- 
tween parameters (see iPal et all 120081) . Times of mid- 
transit for intermediate events were interpolated using 
these two epochs and the corresponding transit number 
of each event, N tr . The eleven main parameters describ- 
ing the physical model were thus T c ._65> T C; o, R p f 7?*, 6 2 , 
K, k = ecosw, h = esinw, K2, faz and hi- Five 
additional parameters were included that have to do with 
the instrumental configuration. These are the HATNet 
blend factor i?i ns t,247, and i?inst,248 which accounts for 
possible dilution of the transit in the 247 and 248 HAT- 
Net light curves from background stars due to the broad 
PSF (20" FWHM), the HATNet out-of-transit magni- 
tudes -Mo, HATNet, 247; and Mo, HATNet, 248, and the relative 
zero-point 7 rc i of the Keck RVs. 

We extended our physical model with an instrumental 



model that describes brightness variations caused by sys- 
tematic errors in the measurements. This was done in a 
simila r fashion to the analysis presented by IBakos et al.l 
(2010). The HATNet photometry has already been EPD- 
and TFA-corrected before the global modeling, so we 
only considered corrections for systematics in the follow- 
up light curves. We chose the "ELTG" method, i.e., EPD 
was performed in "local" mode with EPD coefficients de- 
fined for each night, and TEA was performed in "global" 
mode using the same set of stars and TEA coefficients 
for all nights. The total number of fitted parameters was 
16 (physical model with 5 configuration-related parame- 
ters) + 36 (local EPD) + 10 (global TFA) = 67, i.e. much 
smaller than the number of data points (2866, counting 
only RV measurements and follow-up photometry mea- 
surements) . 

The joint fit was performed as described in lBakos et al.1 
(2010). We minimized x 2 m the space of parame- 
ters using a hybrid algorithm, com bining the down- 
hill simplex method (AMOEBA; see iPress et al.l 1 19921) 
with a classical linear least squares algorithm. Pa- 
rameter uncertainties were derived apply ing the M arkov 
Chain Monte-Carlo method CM CMC. see lFordll2006l) us- 
ing "Hyperplane-CLLS" chains (jBakos et al.ll2010ft . This 
provided the full a posteriori probability distributions of 
all adjusted variables. The a priori distributions of the 
parameters for these chains were chosen to be Gaussian, 
with eigenvalues and eigenvectors derived from the Fisher 
covariance matrix for the best-fit solution. The Fisher 
covariance matrix was calc ulated ana lytically using the 
partial derivatives given bv lPgl(l2009l) . 

Following this procedure we obtained a posteriori dis- 
tributions for all fitted variables, and other quantities of 
interest such as a/R+. As described in Section l3~Tl a/R± 
was used with stellar evolution models to infer a theo- 
retical value for log g* that is significantly more accurate 
than the spectroscopic value. The improved estimate was 
in turn applied to a second iteration of the SME anal- 
ysis, as explained previously, to obtain better estimates 
of Teff* and [Fe/H]. The global modeling was then re- 
peated with updated limb-darkening coefficients based 
on those new spectroscopic determinations. The result- 
ing geometric parameters pertaining to the light curves 
and velocity curves are listed in Table [5] and Table [5] 

Included in Table [5] is the RV "jitter" . This quan- 
tity accounts for RV variability due to rotational mod- 
ulation of stellar surface features, stellar pulsation, un- 
detected planets, and u ncorrected sys tematic errors in 
the velocity reduction (|Wrightl I2005D . Our adopted 
jitter value of 2.3 ms _1 was chosen so that reduced 
X 2 = 1 for the RV data of the global fit. This value 
is consistent with the jitter of an ensemble of chromo- 
spher ically quiet, late G/early K dwarf stars ([Wright 
2005). Auxiliary parameters not listed in the tables 
are: T c _ 65 = 2454449.65967 ± 0.00035 (BJD), T c0 = 
245512L66365 ± 0.00034 (BJD), the blending factors 

Si ns tr,247 = 0.87 ± 0.04 and i?instr,248 — 0.69 ±0.16, and 

7 rc i = 25.0 ± 7.6 ms" 1 . The latter quantity represents 
an arbitrary offset for the Keck RVs, and does not corre- 
spond to the true center of mass velocity of the system, 
which was listed earlier in Table [4] (7rv)- 

The planetary parameters and their uncertainties are 
derived from the a posteriori distributions of the stel- 
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lar, light curve, and RV parameters. We find an in- 
ner planet mass of M p = 0.530 ± 0.018 Mj and a ra- 
dius of R p — 1.010 ± 0.029 Rj, giving a mean density 
p p = 0.64 ± 0.05 gem -3 . These and other planetary pa- 
rameters are listed at the bottom of Table [5] We note 
that the inner planets's eccentricity is significantly non- 
zero: e = 0.346 ± 0.007, lj = 201 ± 2°. 

In addition to HAT-P-17b, we detect a second, outer 
planet in the system. HAT-P-17c is a long-period jo- 
vian planet with a minimum mass ra-i sini2 = 1.4^' \ Mj 
and orbital period Pi = 1798 tgg days. Its eccentricity of 



e 2 



0.101 

o.io * s consistent with a circular orbit. Because 



we have only measured about half of an orbit of HAT-P- 
17c, we conservatively adopt 95.4% confidence intervals 
('2-cr errors') for the error estimates on parameters asso- 
ciated with this planet. (Unless noted, all other param- 
eter uncertainties in this paper are 68.3% confidence in- 
tervals, 'l-cr errors'.) Figure [5] shows the distributions of 
and correlations between m2sini2, P2, and ei from the 
MCMC analysis. Correlations between the Lagrangian 
orbital parameters hi = ei cos W2 and /12 = e.2 sin 0J2 are 
also shown. 

4. DISCUSSION 

We present the detection of HAT-P-17b, a transiting 
hot Saturn in an eccentric orbit, and HAT-P-17c, a cold 
Jupiter near the ice line with an unknown orbital incli- 
nation. In this section we discuss these two planets in 
the context of recent models and trends, the statistics of 
nearly 100 TEPs, and the small number of multi-planet 
systems with one or more transiting members. 

4.1. The Planet RAT-P-llb 

As seen in Figure H3 HAT-P-17b has a radius that is 
typical of other known TEPs with masses in the range 
0.5-0.6 Mj Comparing HAT-P - 1 7b to the theoretical 
models by iFortnev et all ([2007) we find that it is con- 
sistent with gas-dominated planet having a core-mass of 
Mc ~ 25 M e for an age of 4 Gyr, or somewhat less than 
this for older ages. HAT-P-17b is not inflated relative 
to theoretical models. This lack of inflation is consistent 
with the relatively cool temperature (T cq = 792 ± 15 K) 
of HAT-P-17b. 

4.1.1. A Predicted Non-inverted Atmosphere? 

Secondary eclipse measurements of TEPs in multiple 
passbands with the Spitzer Space Telescope have revealed 
two classes of atmospheres among jovian planets. "Non- 
inverted" atmospheres are described well by ID atmo- 
spheric models with water in absorption. In contrast, 
"inverted" atmospheres show dayside emission spectra 
that are best modeled by a high-altitude temperature in- 
version and water in emission. Physically, such an inver- 
sion requires a high-altitude absorber. Planets in inter- 
mediate states (mild inversions , etc) a lso appear possible. 
iKnutson. Howard, k. Isaacson! ((2010) noted that planets 
with inverted atmospheres systematically orbit chromo- 
sphcrically quiet stars while planets with non-inverted 
atmospheres orbit chromospherically active stars. They 
proposed that chromospheric activity, as measured by 
the Ca II H & K index \ogR' HK , traces UV flux which 
is responsible for destroying the photochemically labile 
inversion-causing molecules. These authors are agnostic 



TABLE 5 

Orbital and planetary parameters for HAT-P-17b 

Parameter Value 

Light curve parameters 

P (days) 10.338523 ± 0.000009 

T c (BJD) a 2454801.16945 ± 0.00020 

T w (days) a 0.1691 ± 0.0009 

T12 = T34 (days) a 0.0204 ± 0.0009 

a/R* 22.65 ±0.50 

C/i?* 13.44 ±0.04 

Rp/R* 0.1238 ±0.0009 

0.097i8:°» 

b = acosi/R* 0.31llpog7 

i (deg) 89.21°;'? 

Limb-darkening coefficients b 

ci,i (linear term) 0.3592 

C2,i (quadratic term) 0.2759 

ci,2 0.2843 

c 2 ,z 0.2882 

ci,I 0.3347 

c 2 ,I 0.2800 

RV parameters 

K (ms" 1 ) 58.4 ±0.9 

fc RV c -0.321 ± 0.006 

/i RV c -0.129 ±0.013 

e 0.346 ±0.007 

w (deg) 201 ± 2 

RV jitter (ms -1 ) 2.3 

Secondary eclipse parameters 

T s (BJD) 2454804.246 ± 0.037 

T s 14 0.1326 ±0.0036 

T M 2 0.0154 ±0.0007 

Planetary parameters 

M p (Mj) 0.530 ±0.018 

R p (Rj) 1.010 ±0.029 

C(Mp,R p ) d 0.44 

p p (gem" 3 ) 0.64 ±0.05 

loggp (cgs) 3.11 ±0.02 

a (AU) 0.0882 ± 0.0014 

T cq (K) 792 ± 15 

G 0.108 ±0.004 

F per (10 8 ergs" 1 cm" 2 ) f . . . . 1.95 ± 0.15 

F ap (10 7 ergs- 1 cm- 2 ) f 4.60 ± 0.36 

(F) (10 7 ergs- 1 cm- 2 ) f 8.89 ± 0.66 

a T c : Reference epoch of mid transit that minimizes the correlation 
with the orbital period. It corresponds to Nt r — —32. T14: total 
transit duration, time between first to last contact; T12 — T34: 
ingress/egress time, time between first and second, or third and 
fourth contact. 

k Values for a quadratic law, adopted from the tabulations by 
Claret (2004) according to the spectroscopic (SME) parameters 
listed in Table H] 

c Lagrangian orbital parameters derived from the global modelling, 
and primarily determined by the RV data. 

Correlation coefficient between the planetary mass M p and ra- 
dius R p . 

c The Safronov number is given by — ^ (V^ sc /V or b) 2 — 
(a/H p )(M p /M») (see|H arisen &: Barmar]|2007D . 

Incoming flux per unit surface area, averaged over the orbit. 

on the exact molecule. TiO absor ption has been sup 
gested bv lHubenv et al.l (|2003| ) while lZahnle et all (|2009 
argue that sulfur photochemistry is key to generating an 
atmospheric temperature inversion. HAT-P- 17 is a chro- 
mospherically quiet star with logi?^ K = —5.039. Thus, 
the activity-inversion relation predicts that the atmo- 
sphere ^ofTIATjP^lTbwh^^ The t wo host stars 
in the IKnutson. Howard. &; Isaacson! (|2010D sample with 
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Fig. 5. — A posteriori distributions showing correlations betw 
marked with filled circles. Grayscale regions enclose 68.3%, 95.4%, 

properties most similar to HAT-P-17— WASP-2 (T cff = 
5230 K, \ogR' HK = -5.054) and XO-2 (T eff = 5340 K, 
logi?Jj K = —4.988) — both have planets with measured 
temperature inversions. 

A key caveat is that the activity-inversion relation 
is based on a sample of relatively hot planetary at- 
mospheres. The photochemistry responsible for atmo- 
spheric temperature inversion may not be active on the 
relatively cool HAT-P-17b (T oq = 792 K) For example 
the sulfur photochemistry model by Za hnle et al.l (2009) 
is only active for T > 1200 K. Observations of this planet 
in secondary eclipse will probe the temperature range 
over which such inversions are produced and may offer 
clues to the identity of the absorbing molecules. 

We computed the signal-to-noise ratios (SNR) of warm 
Spitzer secondary eclipse observations of HAT-P-17b 
assuming either efficient day-night circulation (T cq = 
780 K) or a hot day side with no circulation (T cq = 
927 K). For this range of atmospheres, we estimate that 
Spitzer will measure the secondary eclipse depth with 
SNR = 2-5 at 3.6 /im and SNR = 3-6 at 4.5 /im for an 
observation of a single secondary eclipse. This sensitiv- 
ity will help constrain atmosphere models for this planet. 
The superior collecting area and extended IR coverage of 
JWST will give significantly stronger constraints on the 
IR emission spectrum of HAT-P-17b. 

4.1.2. Spin-orbit Alignment 

In the core accretion theory of planet formation, 
hot giant planets like HAT-P-17b form beyond the ice 
line (a few AU from the host star) and subsequently 
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parameters describing HAT-P-17c. Best-fit parameter values are 
99.73% of the MCMC samples. 

migrate inward. Several migration mechanisms have 
been proposed. Tidal intera ctions with the proto- 
planetary disk (jLin et al.lll996| ) deliver gas giants with 

uniformly low obliquit y. Alternatively, Kozai cycles 

(iFabrvckv fc Tremaind I2007D or planet-planet scatter- 
ing (jChattcrie eet al.ll2009D leave the migrated planets 
in high obliquity orbits, possibly with high eccentricity 
(depending on the degree of tidal damping) . Both high 
and low-obliquity systems have been observed by the 
Rossiter-McLaughlin (R-M) effect, suggesting some com- 
bination of migration mechanisms (Fabrycky & Winn 
2 009; Morton fc Joh nson, in preparation). 

I Winn et~ al. (2010) recently noted that nearly all mis- 
aligned (high obliquity) planets orbit hot stars (T e g > 
6250 K). They suggested that all hot giant planets mi- 
grated by one of the high obliquity mechanisms and that 
planets orbiting cool stars subsequently align the spin 
axis of the convective zones and photospheres of their 
hosts with the orbital plane. Stars above this threshold 
temperature lack a significant convective zone and their 
close-in giant planets remain in high-obliquity orbits. 

Although HAT-P-17 is a cool star, the Winn et 
al. model predicts spin-orbit misalignment because the 
wider, eccentric orbit of HAT-P-17b lengthens the 
timescale for orbital decay considerably. Of the known 
TEPs, HAT-P-17b has the longest expected timescale 
except for HD 80606b (J. Winn, personal commu- 
nication). In particular, HAT-P-17b has a longer 
timescal e than WASP-8b, w hich is known to be mis- 
aligned (jQueloz et al.l 120101 ) . This planet (P = 8 d, 
e = 0.31) is broadly similar to HAT-P-17b and also orbits 
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Fig. 6. — Mass-radius diagram of known TEPs (small filled squares). HAT-P-17b is shown as a large filled square. Overlaid are Fortney 
et al. (2007) planetary isochrones interpolated to the solar equivalent semi-major axis of HAT-P-17b for ages of 1.0 Gyr (upper solid lines) 
and 4 Gyr (lower dashed-dotted lines) and core masses of and 25 Mgj (upper and lower lines respectively), as well as isodensity lines for 
0.4, 0.7, 1.0, 1.33, 5.5 and 11.9 gem -3 (dashed lines). Solar system planets are shown with open triangles. 



a relatively cool star (T c ff = 5600 if). 

Measuring the projected spin-orbit angle A of HAT- 
P-17b is a challenging but plausible proposition with 
HIRES. We estimate a 4-11 ms" 1 amplitude R-M effect 
for i>sini = 0.3-0.8 kms -1 . 

4.1.3. Similarity to HAT-P-15b 

HAT-P-17b is strikingly similar to HAT-P-15b 
(jKovacs et al.ll2010l) in orbital period (10.3 d and 10.9 d, 
respectively) and eccentricity (0.35 and 0.19, respec- 
tively). The mass of HAT-P-15b is significantly larger 
though (0.5 Mj and 1.9 Mj, respectively). These two 
planets are the only ground-based transit discoveries 
with orbital periods longer than 10 d. (Other transiting 
planets with P > 10 d include those detected by space- 
based transit surveys and two planets discovered by RVs 
that were later shown to transit.) The bias toward de- 
tecting planets with short orbital periods with ground- 
based transit searches stems from the observational win- 
dow function of a longitudinally-spaced multi-sit e tele- 
scope network (|von Bra un. Kane, fc Ciardil l2009h . The 
HAT- South survey will significantly improve the detec- 
tion of longer period transiting planets wit h a 50% de- 
tectio n rate out to orbital periods of 12 d (jBakos et al.1 
[20091) . 

4.2. Transit Timing Varations 

The presence of a second detected planet in the HAT- 
P-17 system rai ses the possibility of tran sit timing vari- 
ations (TTVs; iHolman fc Murray! [2005h . However, be- 
cause HAT-P-17b and HAT-P-17c are widely separated 



(02/01 ~ 31) and HAT-P-17c is on a nearly circular or- 
bit, the two planets interact very weakly. The TTVs are 
expected to be less than 1 s, undetectable with current 
techniques. 

4.3. The Planet EAT-P-llc 

HAT-P-17c is an approximately Jupiter-mass planet 
separated from its host star by about half the Sun- Jupiter 
separation. Despite having only observed 50% of the or- 
bit of HAT-P-17c, its orbital parameters are well con- 
strained for a long-period planet (Figure [5]). While we 
cannot completely rule out a highly eccentric orbit, only 
3% of the MCMC samples have e > 0.3. 

The 2007 Dec 14 light curve (Figure showing a par- 
tial transit of HAT-P-17b is during the broad transit win- 
dow of HAT-P-17c. We interpret the detected transit as 
due to HAT-P-17b because the timing precisely matches 
the ephemeris derived from other transits. We do not 
detect additional transits (possibly due to HAT-P-17c) 
in that light curve. Because that light curve is the only 
one taken in z-band we cannot compare the photomet- 
ric level of this transit with others to see if it was taken 
entirely when HAT-P-17c was in transit. 

Based on the current orbital fit, the next opportunity 
to search for a transit of HAT-P-17c is in 2012 Oct. 
The timing is favorable for an observing campaign as 
the star is visible for ^5 hr per night from mid-Northern 
latitudes. From the ground, a coordinated, multi-site 
search spanning a range of longitudes is likely necessary 
to rule in or out ~1% deep transits of maximum duration 
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TABLE 6 

Orbital and planetary parameters for HAT-P-17C 



Parameter Value 



RV parameters, as induced by HAT-P-17c 

P2 (days) 1798l|| 

T 2c a (BJD) 2452636 ± 58 

K 2 (ms" 1 ) 25± 2 6 ' 1 

/.-.. o.iH) 

0.02«;j° 9 

^2 o.lo±°;?g 

ui 2 155 ± 150° 

T 2}Peri (days) 2452670 ± 702 

Hypothetical light curve parameters, HAT-P-17c b 

T 2 i4 c (days) 0.883 ± 0.092 

T 2 \i2 = T 34 (days) 0.1008 ± 0.0114 

Hypothetical secondary eclipse parameters for HAT-P-17c a 

T 2b (BJD) 2453531 ± 208 

T 2s 14 (days) 0.92 ±0.11 

T 2a ,i 2 (days) 0.101 ± 0.011 

Planetary parameters for HAT-P-17c 

m 2 sini 2 (Mj) !-4loi 

o 2 (AU) 2.75 ±0.12 

T 2 , eq (K) 140 ± 6 

F 2 , per (lO^rgs" 1 cm" 2 ) d . . 1.1 ±0.5 

F 2 , ap (10 4 erg s" 1 cm" 2 ) d . . . 7.0 ± 1.8 

(F 2 ) (lOVgs" 1 cm" 2 ) d .... 8.7 ±1.2 



a 72c would be the center of transit of HAT-P-17c, if its (unknown) 
inclination is 90°. 

Transits of HAT-P-17c have not been observed. The values are 
for guidance only, and assume zero impact parameter. 
c T14: total transit duration, time between first to last contact, 
assuming zero impact parameter. T12 — T34: ingress/egress time, 
time between first and second, or third and fourth contact. Note 
that these values are hypothetical, and transits of HAT-P-17c have 
not been observed. 

Incoming fiux per unit surface area in periastron, apastron, and 
averaged over the orbit. 



0.883d. 

4.4. Planet Multiplicity 

The migration mechanism of hot jovian planets re- 
mains a major outstanding problem of planet formation 
and evolution. The presence of additional massive plan- 
ets in a system points to migration within the protoplan- 
etary disk, while the absence of additional planets sug- 
gests a more disruptive mechanism such planet-planet 
s cattering or the Koz ai mechanism. 

Wright ct al. (2009) measured the rate of planet mul- 
tiplicity and found that 14% of exoplanet host stars are 
multi-planet systems and another 14% show evidence of 
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multiplicity in the form of an RV trend. Here we com- 
pute the fraction of stars hosting a "cool jovian planet" 
(msini > 0.2 Mj and a > 0.2 AU) that also host a 
"hot jovian planet" (msini > 0.2 Mj and a < 0.2 AU). 
We used the Exoplanet Orbit Databasc(3 of planets with 
well-defined orbital parameters. Of the 375 planet hosts 
(including HAT-P-17), we find 106 stars that host a hot 
jovian planet and 204 stars that host one or more cool jo- 
vian planets. Of the latter group, 10 stars (5%) also host 
a hot jovian planet. Restricting the hot jovian planets to 
a < 0.1 AU, 6/204 = 3% of stars host both cool and hot 

13 http: / /exoplanets.org 

jovian planets. Note that this selection of planets does 
not suffer from a significant detection bias; the Doppler 
signal from a hot jovian planet is essentially always de- 
tectable for systems with a detected cool jovian planet. 
While hot jovian planets represent a disproportionally 
large fraction of the known planets due to observational 
selection effects, multi-planet systems like HAT-P-17 are 
rare. 
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